Chaotically spiking attractors in suspended mirror optical cavities 
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A high-finesse suspended mirror Fabry-Perot cavity is experimentally studied in a regime where 
radiation pressure and photothermal effect are both relevant. The competition between these phe- 
nomena, operating at different time scales, produces unobserved dynamical scenarios where an initial 
Hopf instability is followed by the birth of small-amplitude chaotic attractors which erratically but 
deterministically trigger optical spikes. The observed dynamical regimes are well reproduced by a 
detailed physical model of the system. 
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In the last years, suspended mirror resonators have 
become a research hot topic. In such systems, the intra- 
cavity field and the mirror motion are coupled through 
radiation pressure. In the strong coupling regime they 
could enable the generation of non-classical states of 
light, quantum nondemolition measurements of the field 
quadratures, the creation of entangled states of light and 
mirrors and eventually the observation of the quantum 
ground state of a macroscopic mechanical oscillator (for 
a review of these topics, see Refs. [l], 0]). However, even 
for lower light powers, the action of radiation pressure 
and photothermal effect (i.e., thermal expansion of the 
mirrors due to the absorbed intracavity light) become 
highly nonlinear. In particular, radiation pressure non- 
linearity produces bistability (as shown in the seminal 
experimental works of Refs. 0, Hj]) and, thanks also to 
the delay of the electric field build-up in the cavity, re- 
generative oscillations that eventually become chaotic 0- 
7]. Similar phenomena are generated by thermal effects 
such as photothermal expansion or bolometric force as 
shown, for instance, in Ref. Q. Therefore the interac- 
tion between the two effects, that can both be relevant 
in quantum opto-mechanical experiments, is expected to 
drive the system in regimes of extraordinary complexity, 
that can rise even without the electric field retardation 
@. Indeed, recent theoretical works 1(| U| predict the 
birth of novel kinds of dynamics including spontaneous 
generation of optical spikes. 

Irregular spikes sequences are crucial phenomena in 
the kinetics of chemical reactions 12|, ll3[ and neural dy- 
namics [Til [TBI . While in two-dimensional systems they 
are necessarily related to the presence of noise, in higher 
dimensions they can arise from complicated bifurcation 
cascades, thus having a purely deterministic origin. Typ- 
ical time traces consist of large pulses separated by ir- 
regular time intervals where the system displays small- 
amplitude chaotic oscillations (chaotic spiking). 

In many cases these phenomena can be inter- 

preted in terms of the well known Shilnikov homoclinic 
chaos (l9j . However, Ref. 11 1 proposes that chaotic 



spiking attractors could be observed also in excitable sys- 



tems governed by one slow and at least two fast variables. 
Excitable systems possess a stable attractor (rest state) 
which can be forced to spike by stimuli of amplitude 
above a given threshold |l4|. Each spike is associated 
to a precise orbit in the phase space that do not depend 
on the detailed nature of the stimulus and during which 
the system is insensitive to any new, above-threshold, 
perturbation. If the rest state is a fixed point or a (sub- 
threshold) limit cycle, excitable pulses can only be ob- 
served by adding external perturbations, or in presence 
of noise. On the other hand, a chaotic rest state might 
sporadically overcome the threshold and trigger the ex- 
citable response. The result will be an erratic-though 
fully deterministic-sequence of spikes on top of a chaotic 
background. 

Although similar sequences have been recently repro- 
duced in a semiconductor laser with opto-electronic feed- 
back [2(| , is still lacking an experimental evidence of such 
dynamics in the originally considered system, that in- 
volves peculiar physical characteristics: in particular, the 
multiple timescale competition between dissipative (pho- 
tothermal) and reactive (radiation pressure) nonlineari- 
ties. In this work we report an experimental study of 
such phenomena. 

Our oscillating mirror, realized on a 500 ixm silicon sub- 
strate, consists of a double wheel (see Fig. [T^,b). The 
central mass (main oscillator) and the arms are 70 /im 
thick, while the intermediate ring mass (isolating stage) 
is 500 /.tm thick. On the front side of the wafer, a depo- 
sition of alternate Ta205/Si02 quarter-wave layers pro- 
vides the highly reflective coating. The main oscillator is 
used as end mirror of a L — 12.2 mm long Fabry-Perot 
cavity with a 50 mm radius silica input mirror (trans- 
missivity 7~=110 ppm) operating in a vacuum chamber 
at 5 • 10~ 4 Pa. We obtain a cavity finesse of T = 3 ■ 10 4 
(half linewidth 7 = 225 kHz), with an intracavity optical 
power P c = 1.1 • 10* Pi n , where Pi n is the optical power at 
the input mirror. The mechanical characteristics of the 
micro-mirror has been measured in a previous experiment 
[2lj ]. We find a mechanical frequency wa//27t ~ 250 kHz, 
a mechanical quality factor Q ~ 5 10 3 and an effective 
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FIG. 1: a) Image of the double wheel oscillator, b) FEM 
calculation of the fundamental mode, c) Basic scheme of the 
experimental apparatus. O.I.: optical isolator; H: half-wave 
plate; PD: photodiode; PBS: polarizing beam-splitter; BS: 
beam-splitter. 

mass m ~ 10~ 7 Kg. 

The aim of the experiment is to explore the system 
dynamics as the detuning between laser and cavity res- 
onance, Av, and the injected optical power are varied. 
To this end, the laser is free-running in order to avoid 
possible modifications of the dynamics induced by the 
frequency locking servo-loop. The experimental setup is 
sketched in Fig. [TJ;. The light source is a cw tunable 
Nd:YAG laser operating at A=1064 nm. After a 40 dB 
optical isolator, the beam is injected into the cavity. The 
radiation intensity is controlled by an half- wave plate fol- 
lowed by a polarizing beam splitter and the reflected sig- 
nal beam is monitored by the photodiode PD. 

In Fig. [2] we report the reflected intensity as the laser 
frequency is scanned starting from the red side of the 
cavity resonance. For simplicity, the intensity signals are 
normalized in order to be zero out of resonance and equal 
to 1 at resonance. For low input power we observe the 
typical reflection dip (Fig. For high powers, in- 

stead, the line is strongly asymmetric and dynamical in- 
stabilities arise at both sides of the cavity resonance (Fig. 
[2Jd). In the high-frequency part (blue side), the system 
displays radiation-pressure-induced oscillations with the 
period of the mechanical oscillator, evolving from sinu- 
soidal to chaotic as input power increases. This phe- 
nomenology has been previously observed in toroidal mi- 
crocavities [f|. On the other hand, on the red side of 
the resonance the system exhibits a complex sequence 
of regimes, which have never been observed in such de- 
vices. In particular, we stress that the dynamics is much 
slower than the natural mechanical oscillation and sub- 
stantially different from that observed in Ref. [7[ where 
it is reminiscent of a chaotic damped driven pendulum. 
Figure [2t-g shows five traces of the reflected intensity 
as the laser frequency is delicately increased, approach- 
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FIG. 2: Normalized reflected signal as the laser frequency 
is varied (scanning rate 6.8 MHz/ms) corresponding to a) 
Pin =0.5 mW and b) Pi n ~50 mW. Temporal evolution of the 
reflected intensity for P in =50 mW: c) Av = ~3.9MHz, d) 
Av = -3.64MHz, e) Av = -3.16MHz, f) Av = -3MHz, g) 
Av = — 2.6MHz. Detuning values have been estimated from 
the optical spectrum in Fig. ^jp. 



ing the resonance from the red side. The stationary in- 
tensity value loses stability through a supercritical Hopf 
bifurcation and a small-amplitude quasi-harmonic limit 
cycle is observed (Jit). By further decreasing \Av\, a se- 
ries of bifurcations occurs and the quasi-harmonic limit 
cycle eventually converts into a small amplitude chaotic 
attractor (|2Ji). Although the noise present in the system 
prevents the observation of the full bifurcations sequence, 
the model that we will later introduce suggests that the 
first Hopf bifurcation is followed by a period-doubling 
route to chaos. As the cavity resonance is further ap- 
proached a chaotic spiking regime is observed, where the 
small amplitude chaotic oscillations are sporadically in- 
terrupted by large pulses The mean spikes rate 
increases as |Aj/| decreases (J2J) until the so-called relax- 
ation oscillation regime [3] is finally reached (Hfe). Note 
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that the quasiharmonic cycle and the chaotic oscillations 
are clearly faster than relaxation oscillations. Moreover, 
the large pulses in the chaotic spiking regime have the 
same amplitude and roughly resemble in shape those of 
the relaxation oscillations regime. Indeed, accordingly to 



111 ], each spike should, at least partially, correspond to 
a single relaxation oscillation cycle (i.e., by definition an 
excitable pulse), triggered by sufficiently strong chaotic 
fluctuations. 

The scenario described above is found in a detailed 
physical model of a suspended mirror resonator in the 
presence of radiation pressure and photothermal effects. 
While radiation pressure increases the cavity length by 
pushing the mirrors, the photothermal effect decreases it 
due to their thermal expansion. Therefore, by injecting 
an input field detuned to the red from the cavity reso- 
nance, the first effect tends to increase the intracavity 
optical power in a characteristic time given by 2tt/ljm 
and the second one to reduce it in a time 27r/w c , where 
lu c is the photothermal response angular frequency 22 1. 
For low mechanical quality factors, such competition pro- 
duces excitability and relaxation oscillations in the op- 
tical intensity, through a simple Hopf bifurcation [lOj. 
On the other hand, when the damping rate ujm/Q be- 
comes comparable to lo c , the transitional regime exhibits 
complex bifurcation sequences whose interplay with the 
excitable dynamics leads to a chaotic spiking regime as 



predicted in Ref. [11 1 



At difference with the model analyzed there, in this 
case the optical field evolves on a timescale comparable 
to that of radiation pressure and cannot be adiabatically 
eliminated. The model thus requires to include also the 
intracavity field equation and reads 



-[l-i(S + <p + 9)]e + l , 



as 



9 = -lo ph (6 + 8\s\ 2 ) 



(1) 
(2) 

(3) 



where e is the intracavity field normalized to its reso- 
nant value; <fi and 9 are the cavity length changes due 
to radiation-pressure and photothermal expansion, re- 
spectively, normalized to 7/ = 7—; So is the detun- 
ing between laser and cavity resonance for zero input 
power, normalized to 7; (J — ljm/T and ujph = w c /r are 
the dimensionless oscillator and photothermal frequen- 
cies, normalized to the cavity decay rate T. The pa- 
rameters a — 2 Sj — and B measure the strength of the 
radiation pressure and photothermal [9J effects. Using 
the opto-mechanical parameters given before, we calcu- 
late u) = 1.2 and a = 4.9 10~ 2 Pi n , with Pi n expressed in 
mW. Concerning the photothermal effect, here we remark 
that Eq. ([3]) realistically describes the photothermal dy- 
namics only at frequencies well above ujph- Indeed, the 
frequency response of the length variations to a mod- 
ulated intracavity power has a slow, logarithmic diver- 
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FIG. 3: Numerical integration of Eqs. Q-©, with cj 
Q=5000, a = 4.9 W~ 2 P in , cj ph = 2.4 10 -4 and = 
3.5 10 _1 Pin, with Pin expressed in mW. a) Normalized re- 
flected intensity / = |e| 2 as So is linearly varied in time, 
with a rate of 6.8 MHz/ms, during the numerical integration; 
Pin=0-5 mW (red dashed line) and Pi n — 50 mW (black solid 
line), b) Bifurcation diagram for the peak values of 7 (Pi„=50 
mW). c) Time traces of je| 2 for P in =50 mW: c) S = 4, d) 
S = 4.22, e) S = 4.25, f) <S = 4.29. 



gence at low frequencies [23 , 24 1 , while it follows an intu- 
itive 1 /frequency-dependence as in a standard single-pole 
system at high frequencies. The photothermal response 
is dominated by the fused silica input mirror for which 
we calculate uj c /2it ~ 50 Hz, leading to Dps = 2.4 10~ 4 , 
and 3 = 3.5 10 _1 Pm given an absorption of 2 ppm [9H22]]. 

Fig. [3^, shows the reflected intensity signal as the laser 
frequency is varied. The signal has been obtained by nu- 
merical integration of Eqs. ([TJ-© with Sq linearly vary- 
ing in time. As in the experiment, for high powers the 
line shape becomes strongly asymmetric with dynamical 
instabilities at both sides of the resonance. We remark, 
however, that the line broadening and the width of the 
instability region on the red side are quite different from 
the experiment. This is due to the single-pole description 
of the photothermal effect that, for such slow scanning 
rates, is no longer a good approximation. On the other 
hand, the essential features of the chaotic spiking regime, 
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orbits, associated to a spike in the time series, are always 
approximately equal one to each other. Fig. \4)p shows 
the same reconstruction obtained from numerical data 
showing a good qualitative agreement. 

In conclusion, we have shown that the multiple-time 
scale dynamics induced by radiation pressure and pho- 
tothermal effects in suspended mirror resonators, leads 
to novel dynamical scenarios where a small-amplitude 
chaotic attractor spontaneously generates irregular spike 
sequences. The dynamical mechanism at the basis of 
this phenomenon is relevant also in biology, for instance 
in the study of some apparently erratic neural bursting 
that might have well defined coding functions. Our res- 
onator, operating at cryogenic temperatures and with 
slightly adjusted oscillator parameters, is expected to 
achieve a regime in which quantum fluctuations of ra- 
diation pressure are the dominant noise source (2lj . and 
eventually the oscillator is close to its quantum ground 
state. The instability here investigated, while making the 
observation of quantum-mechanical properties of the sys- 
tem more tricky, might enable the study of the transition 
from classical to quantum domain of a chaotic system 
25]. 

We thank F.S. Cataliotti and M. Inguscio for valuable 
discussion and careful reading of the manuscript. 



FIG. 4: Three-dimensional reconstruction of the phase- 
portrait through the Ruelle-Takens embedding technique, 
from a) experimental and b) numerical time-series of the op- 
tical intensity in the chaotic spiking regime (parameters as in 
Fig. [2^ and Fig. [3b). In both cases we have used a time delay 
t=1.74 us. 



which is significantly faster, are well reproduced. In par- 
ticular, the initial Hopf bifurcation takes place when the 
mean normalized intracavity intensity is approximately 
equal to 0.4 (see Fig. [3b). Then, the system under- 
goes a period doubling cascade leading to the birth of 
small-amplitude chaotic attractors (Fig. [3ji). This is 
illustrated by the bifurcation diagram in Fig. [3}}. Fur- 
ther increasing So, the mean amplitude of the attractors 
grows, until that chaotic fluctuations are sufficiently large 
to eventually trigger excitable spikes (approximately at 
#0=4.23). The result is an erratic - sensitive to initial 
conditions - sequence of pulses on top of a chaotic back- 
ground (Fig. |3k). For larger So, the mean firing rate 
increases until a periodic regime is finally reached. In 
Fig. 2k we plot a three-dimensional phase-space pro- 
jection of the chaotically spiking attractor reconstructed 
from experimental time series by means of the Ruelle- 
Takens embedding technique. This reconstruction shows 
all the essential features of the phase portrait reported in 
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where the small-amplitude chaotic attractor sporad- 
ically triggers large excursions in the phase space. These 
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